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SEQUENTIAL STAR FORMATION IN RCW 34: 
A SPECTROSCOPIC CENSUS OF THE STELLAR CONTENT OF HIGH-MASS STAR-FORMING REGIONS 1 
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ABSTRACT 

In this paper we present VLT/SINFONI integral field spectroscopy of RCW 34 along with 
Spitzer /IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large 
bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class 
II objects is found. At the northern edge of this bubble, an Hn region is located, ionized by 3 OB 
stars, of which the most massive star has spectral type 08.5V. Intermediate mass stars (2-3 M Q ) 
are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. 
North of the Hn region, a photon-dominated region is present, marking the edge of a dense molecular 
cloud traced by H2 emission. Several class 0/1 objects are associated with this cloud, indicating that 
star formation is still taking place. 

The distance to RCW 34 is revised to 2.5 ± 0.2 kpc and an age estimate of 2 ± 1 Myrs is derived 
from the properties of the PMS stars inside the H 11 region. Between the class II sources in the bubble 
and the PMS stars in the H 11 region, no age difference could be detected with the present data. The 
presence of class 0/1 sources in the molecular cloud, however, suggests that the objects inside the 
molecular cloud are significantly younger. 

The most likely scenario for the formation of the three regions is that star formation propagates 
from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, 
after that, the H n region is formed from a dense core at the edge of the molecular cloud, resulting in 
the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular 
cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and 
intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud 
which in turn induces the current star-formation in the molecular cloud, (b) An external triggering 
is responsible for the star-formation propagating from South to North. 

Subject headings: H II regions, stars: formation, stars: early-type, stars: pre-main sequence, photon- 
dominated region (PDR), infrared: stars 
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One of the key questions in astrophysics is under- 
standing the formation and early evolution of high-mass 
stars. They play a major role in shaping the interstellar 
medium due to their strong stellar winds and UV ra- 
diation fields. Observations show that high-mass stars 
preferentially form in clusters; therefore, if we wish to 
understand the formation of high-mass stars, we need 
to understand the clust ered mode of star formation 
(jZinnecker fc Yorkd[2uT51 . 

The actual formation of massive stars happens deeply 
embedded in molecular clouds behind thousands of mag- 
nitudes of visual extinction and is only observable at far- 
infrared to radio wavelengths, providing indirect infor- 
mation about the stars themselves. As soon as massive 
stars begin to emit UV photons and start to develop a 
stellar wind, the surrounding molecular cloud evaporates 
and the stars become visible at shorter wavelengths. 

For the direct detection and the study of the recently 
formed massive stars, the near-infrared is the most suit- 
able regime. The emission of warm dust, the dominant 
source of emission at mid-infrared wavelengths, is much 
fainter in the near-infrared and the actual photospheres 
of the high-mass stars can be probed. The success of this 
approach is demonstrated by several near-infrared imag- 
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ing studies (e.g. iHanson et all [l997HBlum et all 120001: 
iFeldt et al.1 [20031: lAlvarez et alJl2004t IBikl 12004 Kaper 
et al. in prep.). These imaging studies reveal the stel- 
lar content hidden in those star-forming regions, but 
also show the limitations, as these studies are limited 
to regions with Ay less than ~ 30 mag. Near-infrared 
spectroscopy, however, results in a much less ambigu- 
ous identification and classification of the massive star 



content (e.g. 


Watson & Hansonlll997l; iBlum et al. 


2000, 


Hanson et al. 
Puga et al.ll2 


2002; Kendall et al. 2003; Bik et al. 


2005; 


3061). These studies demonstrate the diag- 



nostic power of near-infrared spectroscopy in character- 
izing the young, massive, stellar population. 

Previous near-IR spectroscopic studies only considered 
a limited number of stars per cluster. This incomplete 
spectroscopic sampling of the stellar population, as well 
as the complexity of the fields, prevents any systematic 
investigation of the stellar content in the context of their 
environment. A systematic study on how the stellar con- 
tent depends on the cluster properties like total luminos- 
ity, evolutionary phase as well as their environment re- 
quires a full spectroscopic census of the high-mass stellar 
content. 

We started a program to obtain a full near-infrared 
spectroscopic census of the stellar content of a sample 
of high-mass star-forming regions by means of integral 
field spectroscopy using SINFONI at the VLT. Integral 
field spectroscopy provides a spectrum of every spatial 
element in the field of view. This yields not only a spec- 
trum of all the stars in the observed area, but also emis- 
sion line images of Hil regions, photon dominated re- 
gions (PDRs), and outflows. Supporting Spitzer IRAC 
and MIPS imaging reveal the more deeply embedded ob- 
jects and enable the classification of young stellar objects 
via their colors. 

This paper is the first in a series of studies presenting 
our integral field data set of 10 high-mass star-forming 
regions. As subject of this first study we chose the region 
RCW34, associated with the Vela Molecular Ridge. The 
region consists of a well developed H II region interacting 
with a dense molecular cloud to the North. 

The interaction between the H II region and the molec- 
ular cloud suggests that a triggered star formation sce- 
nario could be applicable to RCW 34. Three different 
mechanisms are discussed in the literature, all result- 
ing from the presence of an expanding Hil region in- 
side a molecular cloud ()Elmegreenll998HDeharveng et all 
2005) . The first mechanism relies on the presence of pre- 
existing molecular clumps. The high pressure of the H n 
region can result in the collapse of those clumps and the 
creation of cometary globules forming a bright rim. The 
second mechanism also produces such globules and bright 
rims, however the clumps are not pre-existing but created 
by sweeping up gas by the H II region. The cores become 
unstable on short time scales and are slowly eroded away 
by the Hn region. The thir d mechanism is the "co l- 
lect and collapse" mechanism (lEhnegreen fc Ladalll977f ). 
Similar to the second mechanism, the H II region sweeps 
up neutral gas; however, in this scenario the neutral layer 
becomes unstable on longer timescales and very massive 
cores are formed. When these massive cores collapse they 
form new star clusters. 

We will determine the properties (spectral-type, mass, 
age) of the high- and intermediate mass stellar popula- 



tion of RCW34 by means of near-infrared integral field 
spectroscopy and Spitzer /IRAC images. The cluster en- 
vironment is characterized via Hn, He I, and H2 emission 
lines. The three triggering mechanisms are tested on the 
observed properties of RCW 34 and it turns out that an 
unusual mode of triggering might be at work here. 

1.1. RCW 34 

The Hn region RC W 34 (IRodeers et al.lfl960h (IRAS 



08546-4254 or Gum 29 (!Gjim|l955])) is located in the Vel a 



Molecular Ridge (VMR) cloud C (|Liseau et al.1 11992T) . 
The region is an optically visible H 11 region of about 2' x 
2' in extent and slightly offset from the main dense cores 
defining the VMR cloud C detected by CO ob servations 
(jMurphv fc Mavl 119911; lYamaguchi et all 119991) . Due to 
a similar velocity (vj gr = 6 km s _1 ), iMurphv fc Mavl 
(|1991l ) argue that RCW34 is likely part of the C-cloud. A 
similar velocity has been measured from the CS line (vj sr 
= 5.5 km s" 1 . Bronfma n et al.l[l9 96D. confirming the as- 
sociation with cloud C, which would place the region at 
a distance somewher e between 1 and 2 kpc. 

Herbst (1975) and Ivan den Bergh fc Herbstl (fl97l al- 
lowed that the exciting star of RCW 34 might be more 
distant, but favored an association with the Vela-R2 
R-association at a distance of 870 pc. Optical spec- 
troscopy by Heydari-Ma layeri (1988) implies that the 
central source vdBH 25a (Ivan den Bergh fc Herbsd[l975l) 
is of spectral type 08.5V. This suggests a distance of 2.9 
kpc, supporting the idea RCW 34 is more distant than 
Vela-R2. In §4.1 we use near-infrared spectroscopy of the 
cluster OB stars to derive a distance in rough agreement 
with Heydari-Malayeri (1988). Heydari-Malayeri also de- 
tected a North-South velocity gradient in the ionized gas, 
which he suggested to be evidence that RCW 34 is an HII 
region undergoing a champagne flow. This is supported 
by Deharveng et al. (2005) who conclude, based on MSX 
and near-infrared observations, that the region is not the 
result of the collect and collapse triggered star formation 
scenario. 

Several signs of active star formation are present 
in RCW 34. Two water masers have been reported 



(|Braz & Epchteinlll983UCaswell et al. 


1989J) as well as a 


class II methanol maser (Chan et al.l 


1996: Walsh et al.l 



119971 ). Associa ted with the water m asers, CS emission 
has been found dZinchenko et al.lll995l) with an estim ated 
mass of 1320 M ? . iGvulbudaghian fc Mavl (|2007t) re- 
ported the detection of 1.2 mm continuum emission to- 
wards RCW34. 

The outline of the paper is as follows; Section 2 outlines 
the data reduction of the SINFONI and Spitzer /IRAC 
data sets used in this paper. In section 3, the interaction 
between the H 11 region and the cluster surroundings is 
studied using the IRAC data as well as the SINFONI 
emission line maps. Section 4 focuses on the stellar spec- 
tra obtained in the SINFONI field; spectral types, lu- 
minosity class as well as extinctions are derived, and in 
Section 5 the stellar content is discussed and an age es- 
timate is obtained for the NIR cluster in RCW 34. In 
Section 6 we discuss the complete region and combine 
all the different properties to derive a possible formation 
scenario for RCW 34, and Section 7 summarizes the main 
conclusions of the paper. 

2. OBSERVATIONS AND DATA REDUCTION 
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2.1. SINFONI observations 

Near-infrared H- and K-hsaid observations were per- 
formed using the Integ r al Field instrument SINFONI 
(jEisenhauer et all 12001 iBonnet et all [2004) on UT4 
(Yepun) of the VLT at ESO, Paranal, Chile. The ob- 
servations of RCW 34 were performed in service mode 
between February 4 and 9, 2007. The non-AO mode of 
SINFONI, used in combination with the setting provid- 
ing the widest field of view (8"x8"), delivers a spatial 
resolution of 0.25" per slitlet. The typical seeing during 
the observations was 0.7" in K. The H+K grating was 
used, covering these bands with a resolution of R=1500 
in a single exposure. 

We observed RCW 34 with a detector integration time 
(DIT) of 30 s per pointing to guarantee a good S/N (~ 
70) for the early B-type stars within this cluster. RCW 
34 is much larger than the field of view of SINFONI, 
the near-infrared extent of the Hn region is 104" x 60" 
on the sky. The observations were centered on the cen- 
tral 08.5 V star (vdBH 25a) at coordinates: a(2000) = 
08 h 56 m 28 s . 1, (5(2000) = -43°05'55.6". We observed this 
area with SINFONI using a raster pattern, covering ev- 
ery location in the cluster at least twice. The offset in 
the East- West direction was 4" (i.e. half the FOV of the 
detector) and the offset in the North-South direction was 
6.75". A sky frame was taken every 3 minutes using the 
same DIT as the science observations. The sky positions 
were chosen based on existing NIR imaging in order to 
avoid contamination. Immediately after every science 
observation, a telluric standard of B spectral type was 
observed, matching as close as possible the airmass of 
the object. This star is used for the removal of the tel- 
luric absorption lines as well as for flux calibration. The 
stars chosen as telluric standards were HIP 048652, HIP 
045117, HIP 052202 and HIP 040629, one for each of the 
observing blocks. 

2.1.1. Data reduction 

The SINFONI data were reduced using the SPRED 
(version 1.37 ) software developed b y the MPE SINFONI 
consortium (jSchreiber et all l2004t lAbuter et aTll2006h . 
This software was used to apply all the basic steps of 
the data reduction like flat-field correction, dark removal, 
distortion correction, wavelength calibration and correc- 
tion of bad pixels. Special care was taken to remove bad 
pixels from the calibration frames. The construction of 
the final calibration data was done in three iterations by 
first constructing a preliminary calibration set used to 
detect and remove the bad pixels using the 3D informa- 
tion. The cleaned raw files are used to construct a new 
calibration set and the procedure is repeated. 

The bad pixels on the science data are removed us- 
ing the same 3D interpolation method. However, due to 
severe cross talk between the pixels, the surrounding pix- 
els are affected at lower intensity and are not properly 
removed by the algorithm. Therefore we apply an addi- 
tional bad pixel removal later in the reduction process. 
The cleaned science and sky frames are reconstructed 
into 3D cubes using the calibration data. Also the dark 
frames are subtracted and the cubes are corrected for 
atmospher ic diffraction. We applied the procedure de- 
scribed by iDaviesI (|2007t ) to remove OH line residuals. 

In the reduced cubes we discovered a drift of the stellar 



position versus wavelength. The shift is strongest in K 
(up to 2 pixels) and resulted in an elongation of the stars 
in the collapsed image. We measured the position of 
the telluric standard star versus wavelength using a 2D 
gaussian and fitted a spline to the change of position in 
both spatial directions. This spline fit was used to correct 
the science cubes. The residual measured on some bright 
stars in the science frames was within 0.5 pixels, which 
is acceptable for our purposes. 

The last step in the reduction process is flux calibra- 
tion and removal of the telluric absorption lines. This 
is done using the extracted spectrum of a B type stan- 
dard star, using an aperture radius of 6 pixels (0.75"). 
After the extraction, the correction for the telluric ab- 
sorption lines in the science spectra is done in two steps: 
First, the telluric standard star is corrected using a high 
signal-to- noise atmospheric template spectrum 15 . This 
high S /N spectrum is not taken at the same airmass and 
atmospheric conditions as our science spectra, but pro- 
vides a first rudimentary correction in order to reduce the 
absorption contamination of the telluric standard star. 

Second, the Brackett and He I lines are fitted with a 
Lorentzian profile. The line-fits are subtracted from the 
original telluric standard star spectrum, leaving an at- 
mosphere transmission spectr um. The flux calibr ation of 
the spectra uses the 2MASS (|Cohen et al.ll2003D magni- 
tude of the standard stars. As we only use normalized 
spectra for the point sources and flux ratios for the emis- 
sion line spectra, the uncertainty in the flux calibration 
due to atmospheric variations is not critical. 

The correction of the atmospheric absorption lines and 
the flux calibration are performed on the 3D cubes; the 
spectrum of every spatial pixel is divided by the spectrum 
of the telluric standard star and flux calibrated. Finally, 
the individual 3D cubes of one observing block (between 
50 and 60 cubes) were combined into one large 3D cube. 
For RCW 34, four of those large cubes are needed to 
cover the near IR emission of the cluster. 

Before extracting the point source spectra from the sci- 
ence cubes, an attempt was made to remove the extended 
background emission from the data. This background 
emission is primarily nebular emission from the HII re- 
gion and the PDR (Brackett lines, Hei, H2). A median 
filter with a kernel of 5 pixels is used to remove the ex- 
tended nebular component in the cubes. We applied this 
median filter 3 times to remove any high frequency vari- 
ations present which could be due to stars. This median 
filtering has only a small effect on the stellar continuum; 
the contribution in the total flux is less than 1 %. The 
point sources were extracted using a circular aperture 
with a radius of 6 pixels, analogous to that of the stan- 
dard star. 

Several maps centered on emission lines as well as con- 
tinuum maps are created from the combined SINFONI 
data cubes. The maps are made by summing 4 to 7 
slices of the data cube in the wavelength direction. The 
images in the spectral lines are corrected for continuum 
emission from the stars by interpolating each pixel of the 
adjacent continuum images (bluewards and redwards of 
the line, respectively). This results in a clean subtrac- 
tion for most of the stars. However, for the brightest 

15 taken from the SINFONI website and obtained at NSO/Kitt 
Peak Observatory 






Fig. 1. — Spitzer /IRAC 3 color image of RCW 34 and its surroundings (blue: 3.6 (im, green: 4.5 (im, red: 8 (im, color version online). 
The top left corner is only covered at 3.6 fim. The area delimited by the white box is the area covered by the SINFONI observations. The 
location of the central star of RCW 34, vdBH 25a (a(2000)= 08 h 56 m 28 s . 1, (5(2000) = -43°05'55.6".), is shown with an asterisk. The Hn 
region surrounds vdBH 25a and is fully covered by our SINFONI observations. 



stars in the maps, residuals remain present. An interpo- 
lation of the line emission a bit further away from the 
stars successfully restored these residuals. 

We created velocity maps of the strongest emission 
lines, Br7 and H2, for every (spatial) pixel in the data 
cube. The central velocity of the Br7 and H2 emission 
line is derived from a gaussian fit to the line profile. To 
obtain better accuracy in the wavelength calibration we 
performed the same fitting procedure to a nearby OH 
skyline. Finally, the velocities were converted to the local 
standard of rest. The errors were calculated by quadrat- 
ically summing the fitting errors on the nebular line and 
the OH line and are on the order of 10 - 15 km s _1 . 

2.2. IRAC data reduction and photometry 

Additional data obtained with IRAC (Fazi cTet al.l 
120041) on board of the Spitzer satellite have been retrieved 
from the Spitzer archive (PID: 40791, PI: Majewski). 
The data were obtained on January 31, 2008 in a map- 
ping mode, covering a large part of the outer galaxy. We 
selected the frames covering the area around the cluster 
as well as an adjacent control field. Each frame has an 
exposure time of 1.2 seconds and every position is cov- 
ered twice. The raw frames were processed using the 



standard pipeline version S17.02 to produce the basic 
calibration data set. This dat a set was processed fu rther 
using MOPEX version 16.2.1 (jMakovoz et al.l 120061 ). 

As the IRAC data are undersampled, aperture pho- 
tometry is preferred above psf-fitting photometry. Source 
detection and photometry were performed using the 
Daophot package under IRAF. Point sources with a peak 
level of more than 5 times the RMS of the background 
were detected. Aperture photometry was applied with an 
aperture of 3 pixels (3.7") and a sky annulus between 3 
(3.7") and 7 (8.5") pixels. These small values were chosen 
due to the variable extended emission. The appropriate 
aperture corrections and zeropoints were adopted (Ta- 
bles 5.1 and 5.7 of the IRAC Data Handbook) to obtain 
the magnitudes of the objects. 

3. CLUSTER SURROUNDINGS 

NIR imaging (Bik 2004; Deharveng et al. 2005) shows 
an H 11 region extending about 1.5' x 1' around the central 
star vdBH 25a. The IRAC images (Fig. 1) show a much 
larger star-forming region than is visible in the NIR. The 
region extends to the South and East and shows an open 
bubble structure. The total size of the mid-IR region is 
about 6' x 9'. This structure looks similar to the (partial) 
bubbles identified in the GLIMPSE survey (Churchwell 
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0.5 1.0 
[5.8]-[8.0] 

Fig. 2. — Spitzer /IRAC color-color diagram. The star-shaped 
points represent stars without infrared exc ess, class II stars ar e 
plotted as triangles inside the squared box (Mcgcath ct al. 2004). 
The diamond symbols represent the location of the protostars (class 
0/1) and the squares show the location of the class I/II sources. 
The numbered objects are located inside the SINFONI field and 
are discussed in more detail in Sect. 4.3. 

et al. 2006, 2007). 

In this section we will determine the nature of the stars 
located inside the bubble using the IRAC photometry 
and investigate whether the objects in the bubble are re- 
lated to those detected in the H II region. We analyze the 
extended emission seen in the SINFONI observations and 
derive properties of the H II region and the surrounding 
molecular cloud. 

3.1. Spitzer/IRAC imaging of the cluster and its 
surroundings 

A 3-color Spitzer /IRAC image of the region is shown 
in Fig. [Hand Fig. [2] shows a [3.6]-[4.5] vs. [5.8]-[8.0] 
color-color diagram of the stars in the field. Objects not 
detected in all bands are not included in the analysis. 
We have not plotted the stars in the control field which 
only contains star s without IR excess. W ith t he color cri- 
teria p roposed by iMegeath et al.1 ()2004h and lAllen et aD 
(2004) we can identify several classes of objects. Stars 
without IR-excess have colors close to 0.0. The IRAC col- 
ors are insensitive to the effective temperature of the pho- 
tospheres. These objects, plotted as star-shaped sym- 
bols, could be fore- or background late-type stars as well 
as earlier-type main-sequence stars associated with the 
star-forming region. The objects with [3.6]- [4.5] between 
0.0 and 0.8 and [5.8]-[8.0] between 0.4 and 1.1 are clas- 
sified as class II objects. This region is marked with a 
box, and the objects showing class II colors are plotted as 
triangles. Objects with [3.6]-[4.5] > 0.8 and [5.8]-[8.0] > 
0.0 or [3.6]-[4.5] > 0.4 and [5.8]-[8.0] > 1.1 are classified 
as protostars (class and I objects, diamonds), objects 
with [5.8]-[8.0] > 1.1 and [3.6]-[4.5] between 0.0 and 0.4 
are classified as class I/II sources (square symbols). The 
5 numbered objects are located in the SINFONI field and 
are discussed in more detail in Sect. 4.3. and 5.3. 

The cl ass I/II identification s hould be considered with 
caution. IMegeath et al.1 (|2004f ) show that these objects 
are likely class II sources with the 8 /im flux contam- 
inated by variable extended emission. Fig. Q] shows 
that these objects are located on top of a strong dif- 
fuse background. Therefore, we consider them as class II 
objects. Additional sources of contamination, especially 



for the detection of protostars, could be planetary neb- 
ula, A GB stars or background galaxies (IWhitnev et al.1 
2003). However, the absence of sources with these col- 
ors in the control field suggests that the contamination 
is negligible. 

We have plotted the different classes on the IRAC 3.6 
/im image (Fig. [3]) to study their potential relation with 
the extended MIR emission of the bubble. The locations 
of the stars without infrared excess (white, star-shapes) 
are not related to the extended emission, but evenly dis- 
tributed over the field outside the bubble. These objects 
are most likely field stars. 

The location of the YSOs shows a clear correlation with 
the bubble and the extended emission. Most of the class 
II sources (yellow triangles) are located inside the bubble. 
The protostars (red triangles) are spatially associated 
with the extended MIR emission. Most of these objects 
are located near the bright arc in the north, but also some 
objects are associated with the fainter extended emission 
in the southern part of the bubble. Although no obvious 
cluster is visible inside the bubble, the class II objects 
located inside the bubble suggest a recent star formation 
event. This mi ght have been responsib le for the creation 
of the bubble (|Churchwell et all 120061 ). The location of 
the younger protostars in the outer regions of the bub- 
ble indicates the presence of a more recent generation of 
stars near the edge of the bubble. 

3.2. SINFONI line and continuum maps 

The SINFONI observations are centered on the ioniz- 
ing source, vdBH 25a, and cover both the H n region as 
well as the bright arc seen in the IRAC images. Fig. 2] 
shows the 3 color composite of 3 linemaps along with con- 
tours of the IRAC 3.6 /im image. This image reveals the 
stellar content of this region and the overlay shows that 
several near-infrared sources are also detected at IRAC 
wavelengths. The region is dominated by the two bright 
stars in the center of the field (one of them is vdBH 25a 
at the location asterisk symbol) and a relatively sparse 
cluster of about 130 stars down to K=17 mag, surround- 
ing the central region. For all the point sources visible in 
this image a SINFONI H- and .ftT-band spectrum is avail- 
able and a spectral classification is obtained (Sect. 4). 
The Hn region is traced by Br7 (2.16 /im, red) and Hci 
(2.058 /im, blue) emission and the surrounding PDR by 
the molecular H2 emission (2.12 /im, green). Just North 
of the position of the central star, the Br7 emission is 
strongest. Even further to the north, the Br7 and Hci 
emission drop abruptly and the H2 emission begins, in- 
dicating the interface between the ionized and molecular 
gas, coincident with the bright arc seen in the IRAC im- 
ages. 

3.2.1. H II region 

We have extracted spectra at six locations covering 
most of the Hn region (annotated with numbers I - VI 
in in Fig. 2J. Table [T] gives the position and size of the 
six extracted regions. 

As an example, the upper panels of Fig. [5] show the 
spectra extracted from region II, just North of the ioniz- 
ing star. 

To check the hypothesis of a champagne flow, we esti- 
mate the column density of the molecular cloud in which 
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Fig. 3. — Spitzer /IRAC 3.6 /an image of RCW 34 and its surroundings. The different classes of objects as derived from the IRAC 
color-color diagram are overplotted. The symbols used are identical to those in Fig. [2] 




32.0 31.0 8:56:30.0 29.0 28.0 27.0 26.0 25.0 24.0 

Right ascension 

Fig. 4. — Three-color image created from the SINFONI line maps of RCW 34 (color version online). Red: Br7 (line + continuum) 
emission, green: H2 (2.12 fim, line + continuum) and blue: He I (2.058 fim, line + continuum). The horizontal striping is a data reduction 
effect. Overplotted in contours is the 3.6 (im image. The letters b, c mark the parts of the H2 emission (region (a) is not annotated but 
comprises the entire arc North of the H II region) for which excitation diagrams are constructed to identify the nature of the emission (Fig. 
[BJ. The six boxes show the regions from which spectra are extracted to derive the extinction and the Hel/Br7 ratio. The asterisk shows 
the location of the main ionizing source vdBH 25a. 
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the H II region is embedded. In an H n region, the rela- 
tive line intensities of the hydrogen recombination lines 
are mostly determined by the Einstein coefficients, with 
a very weak dependence on dens ity and electron temper- 
ature (jStorev fe Hummerl [1995) . For example, the ra- 
tio of the BrlO (1.73 /im) and the Br7 (2.16 /im) lines 
is predicted to be 0.33. Deviations from this value are 
caused by foreground extinction. We have converted the 
measu red Brl0/Br7 ratio to Ay using the lCardelli et al.l 
( 1989) extinction law. The derived Ay values for the six 
extracted regions are listed in column 7 of Table [T] The 
bulk of the H II region has an Ay of around 5-6 mag (re- 
gions I, II and IV). In region III, on the border with the 
molecular cloud, we measure an Ay of 15.9 mag, while 
towards the South the extinction declines to 4.4 mag and 
1.5 mag for regions V and VI respectively. This extinc- 
tion gradient can be explained as a density gradient in 
the surrounding molecular cloud and is consistent with 
the presence of a champagne flow. 

The (extinction corrected) ratio of the 1.70 /im He I line 
over Br7 is dependent on the temperature of the ioniz- 
ing star and can be used to constrain its spectral type 
(jLumsden et al. 2003|). The extinction correction line ra- 
tios for the six regions are given in column 6 of Table [TJ 
The maximum ratio is measured just South of the loca- 
tion of the ionizing stars, and has a value of 0.047 ± 0.004 
in region V and 0.06 ± 0.01 in region VI. Compar ing this 
value to the predictions of ILumsden et all ()2003D . an ef- 
fective temperature for the central star around 35000 K is 
expected, which corresponds to a sp ectral type between 
08 V and 09 V (jMartins et al.ll2002l ). c onsistent with the 
optica l spectral type determination of IHevdari-Malaveril 
(fl988l) . 

3.2.2. PDR and outflows 

The IRAC detections of class and I sources show that 
stars are still forming close to the bright arc. An addi- 
tional argument for active star formation is the presence 
of outflows, detectable by H2 line emission (bottom pan- 
els Fig. [5]). 

The H2 emission in the near-IR can be either ther- 
mally or non-thermally excited. The thermal emission 
originates from shocked neutral gas being heated up to a 
few 1000 K (outflows), while the non-thermal emission is 
due to fluorescence by non-ionizing UV photons between 
912 and 1108 A. These two mechanisms can be distin- 
guished as they preferentially populate different energy 
levels giving rise to different line ratios. The fluorescence 
mechanism preferably excites the higher vibrational lev- 
els, wh ile for thermal exc itation the lower levels are pre- 
ferred (jDavis et al.ll2003ft . 

Fig. 2] shows the location of the 2.12 /im H2 emission 
(green) in relation to the emission from ionized gas (B1-7 
and He 1). The H2 emission has a very filamentary struc- 
ture and appears just North of the ionized gas, co-located 
with the bright arc in the IRAC images. Additionally, 
two compact sources are visible in the north-east. To 
identify the excitation mechanism of the arc of filamen- 
tary emission (region (a)) as well as that of the two com- 
pact regions ((b) and (c), top left panel of Fig. [6]) we 
construct excitation diagrams for these regions. 

First, the line flux ratios are extinction corrected. The 
ratio of the 1-0 S(l) (2.12 /im) over the 1-0 Q(3) (2.423 
/im) is insensitive to the ionization mechanism and has 



an expected value of 0.7 ((Turner et al.lll977l ). Therefore, 
the observed ratios are used for an extinction estimate in 
the same manner as discussed in Sect. 3.2.1 (Table [2]). 
Region (a) covers a large area, 262 square arcsecs and 
shows a large extinction variation. The derived prop- 
erties are dominated by the strongest emission regions, 
located close to the Hil region where the extinction is 
lowest. 

Fig. [6] also show the excitation diagrams of the 3 re- 
gions. In these diagrams, the measured column densities 
of the lin es are plotted against the energ y of the upper 
level (see iMartin-Hernandez et al.l ((2008) for a detailed 
description of this diagram). The total column den- 
sity w as calculated using the description of iZhang et al.l 
( 1998). Different vibrational levels are indicated with dif- 
ferent symbols. Region (a) has the most lines detected 
as it covers the largest area (the complete arc). For the 
spectra of regions (b) and (c) , the weaker lines seen in the 
spectrum of region (a) are not detected and 3 a upper 
limits are given instead. The solid line in the diagrams 
is a single temperature fit to all the data points, while 
the dashed line in the diagram of region (a) represents 
a linear fit to only the 1-0 S lines. For region (a) it is 
clear that the column densities are not represented by 
a single temperature gas, but that the line fluxes of the 
2-1 and 3-2 vibrational levels are higher than expected 
from the 1-0 line fluxes. A likely excitation mechanism 
of the H2 gas in the arc is fluorescence by non-ionizing 
UV photons. We have searched for spatial variations of 
the temperature and column density in region (a) by cal- 
culating the excitation diagrams for several sub-regions, 
but could not detect any significant deviations from the 
values of the total area. 

The excitation diagrams of regions (b) and (c) are well- 
represented by a single temperature. We note, however, 
that the 3-2 lines are not detected in either region, and 
the 2-1 lines are mostly upper limits. We conclude that 
the excitation diagrams of both regions are consistent 
with shock excited emission in an outflow. Their elon- 
gated shapes, especially for region (b), are also suggestive 
of outflow. The linear fits of the excitation diagrams also 
allow us to determine the column density and tempera- 
ture of the emitting gas (Table [5]) . 

An additional constraint on the nature of the H2 emis- 
sion regions is the velocity of the different regions. The 
obtained velocities are listed in the last column of Table 
[21 The two outflows have a clear blue shifted velocity 
compared to that of the extended PDR, with region c 
having the largest value of — 26 ± 16 km s _1 . These rel- 
atively low velocities as well as the absence of [Fen] at 
1.644 /xm suggest tha t both shocks are C-type shocks 
(jStahler fe Pallall2005h . The more destructive shocks (J- 
type) would hav e strong [Fe III emission as w ell as larger 
shock velocities (jHollenbach fe McKedll989f ). The veloc- 
ity measured in the arc (region (a)) of 1 ± 14 km s _1 
is more blue shifted than the velocity of the H 11 region, 
but consistent with that of the Vela molecular cloud: 6 
km s -1 . 

For outflow (b) a driving source could not be identified 
in the IRAC images, most likely the source is deeply em- 
bedded and not visible at 8 /im. For source (c), however, 
we could identify a possible driving source. About 3" 
South-East of the H2 emission a deeply embedded proto- 
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TABLE 1 
Properties of the Hii gas 



Region 


R.A. (2000) 
(h m s) 


DEC. (2000) 
(° ' ") 


size 
(") 




(km s _1 ) 


Hel/Br7 


A v 
(mag) 


I 


08 56 29.8 


-43 05 52.4 


2.4 X 


4.6 


18 ± 13 


< 0.13 


5.3 ± 0.2 


11 


08 56 28.1 


-43 05 50.8 


13.7 x 


1.9 


16 ± 8 


0.038± 0.002 


5.1 ± 0.1 


111 


08 56 27.3 


-43 05 39.2 


2.3 x 


2.7 


16 ± 9 


< 0.17 


15.9 ± 0.6 


IV 


08 56 25.7 


-43 05 55.7 


5.0 x 


3.0 


10 ±12 


0.032 ± 0.006 


5.6 ± 0.2 


V 


08 56 27.7 


-43 06 05.0 


7.7 x 


4.0 


13 ± 9 


0.047 ± 0.004 


4.4 ± 0.1 


VI 


08 56 27.7 


-43 06 11.9 


7.7 x 


3.7 


11 ± 10 


0.06 ± 0.01 


1.4 ± 0.2 
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TABLE 2 
Physical properties of the H2 gas. 



Region 


Trot 


N(H 2 ) 








(K) 


(cm" 2 ) 




(mag) (km/s) 


a 


3720 ± 90 


3.6 ± 0.6 x 


10 17 


9 +1 ± 14 


b 


2217 ± 90 


1.8 ± 0.4 x 


10 18 


13 -14 ± 14 


c 


1980 ± 115 


4.6 ± 2.0 x 


10 18 


20 -26 ± 16 



Fig. 5. — Top panels: H- and A"-band SINFONI spectra of the Hii region in RCW 34. The spectra are drawn from a region just North 
of the ionizing stars (region II in Fig. [4). The spectra show numerous hydrogen lines of which Br7 (2.166 (im) is the most prominent. Also 
several Hel lines are visible, suggesting that the central star is of spectral type late O. Bottom panels: Spectra of the photo dominated 
region (PDR), North of the Hii region (region (a)). The spectrum is dominated by H2 lines. Also some faint [Fe II] emission is present. 
This line can be formed in J-shocks and also in the low density material between the H II region and the PDR. 

vdBH 25a (08.5 V). The detected North-South extinc- 
tio n gradient as well as the Ha velocity gradient found 
bv iHevdari-Malaveril (|1988h suggest that the H II region 
is a champagne flow. 

In between the molecular cloud and the H 11 region a 
PDR is located traced by H2 emission and 8 fim emission. 
The star formation inside the molecular cloud seems to 
be more recent as masers, class 0/1 sources and outflows 
are present in that area. 

4. SINFONI NEAR-INFRARED SPECTROSCOPY 

This section focuses on the stellar content of the NIR 
cluster inside the Hn region RCW 34. In total 134 stars 
have been identified between K—8.7 and 17 mag (see 
Fig. 2] for their location). Our SINFONI observations 
provide an H- and i-T-band spectrum of each source with 
a spectral resolution of R=1500. Stars with K < 14 have 
high enough S/N spectra to obtain a reliable spectral 
type. This reduces the sample to 26 stars, which will be 
discussed in the remainder of this paper. 

To estimate the stellar content, Fig. [7] shows the K vs. 
(J-K) color-magnitude diagram (CMD) of the area i n 
RCW 34 covered with SINFONI taken fromES IpOOl . 
who imaged a 5'x 5' field using SOFI at the NTT tele- 



star is identified on the IRAC images with a very bright 
4.5 /xm flux (source # 28 in Fig. [5]). The location of 
the source is shown in Figs. [4] and [6] where the 3.6 /zm 
contours of source # 28 are overlaid on the SINFONI 
images. 

The excess 4.5 /im emission is most likely coming from 
H 2 emission closer to the source, too deeply embedded 
to detect at near-infrared wavelengths. 

3.3. Summary 

Summarizing, RCW 34 consists of 3 distinct regions. 
A large bubble is detected in the IRAC images, with 
a cluster of class II stars inside. At the northern edge 
of this bubble the Hn region is located, ionized by 
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Fig. 6. — Top left: zoom on the SINFONI H2 map highlighting the two candidate outflow regions (b) and (c). The contour is the 3.6 (im 
image highlighting source #28, the possible driving source for region (c). Panels 1 (region (a)), 2 (region (b)) and 3 (region (c))show the 
excitation diagrams of 3 different regions in the PDR identified in Fig [4] The solid line in all the diagrams is a single temperature fit to all 
the lines detected, while for the dash-dotted line in panel a only the 1-0 S lines are included in the fit. The excitation diagram of region 
(a) is that of the complete arc North of the H II region. The excitation diagram suggests that UV fluorescence is the excitation mechanism 
in this area. Regions (b) and (c) are most likely outflows as their excitation diagrams suggest thermal excitation. 



scope. 

The symbols in the plot represent the different spec- 
tral types derived later in this section. The CMD shows 
that the cluster consists of stars showing a large range in 
J — K colors, which could be related to large extinction 
variations in the cluster. In the remainder of the paper, 
the brightest star (vdBH 25a) is labeled #1 and the la- 
bel number increases with if-band magnitude. Table |5] 
provides the list of the 26 brightest stars, as well as some 
additional objects of special interest, and their properties 
derived in this paper. 

4.1. Early type stars 

The SINFONI H+K spectra of the three brightest clus- 
ter members (sources #1, #2 and #3) show Br7 absorp- 
tion in the if-band and BrlO - 14 absorption in the H- 
band (Fig. [3]). Weaker lines of He I are also visible in 
H (1.70 fim) and K (2.058 and 2.11 /im). To obtain 
a spectral type estimate we appl y the if-band classifi- 
cation scheme for OB stars from lHanson et all (jl996). 
Their classification links a K-band spectral type to an 
optical spectral type. For late O and early B stars, the 
classification is based on Hei (2.11 /im) and Br7. For 
the i? -ban d, we used the cla ssification of lHanson et al.l 
(1998) and lBlum et al.l (|1997l) where the important lines 
are Hei (1.700 //m) and Brll (1.68 fim). 



Table [3] shows the measured Equivalent Widths (EW) 
of the relevant lines in the H- and if -band spectra of 
the 3 OB stars. In column 6 the if-band spectral type is 
given, with the corresponding optical spectral type listed 
in column 7. The brightest star, #1, has a spectral type 
between 08V and B1V based on the presence of Hei 
(2.11 /xm) and moderately strong Br7. Star #2 has a 
much stronger Br7 absorption and is therefore classified 
as B1V - B2.5V. The fainter OB star (#3) does not have 
Hei (2.11 fim) in its spectrum and is therefore classified 
as B3V - B7V. The EWs of the iJ-band lines are in agree- 
ment with the if-band spectral type. The only deviating 
value is the very strong Brll absorption of #3, which is 
about equally strong as the Br7 line. The strength of 
Brll is compatible with that of a supergiant. 

The spectral type of star #1 (08V-B1V) is compat- 
ible with the opti c al cla ssification as an 08.5V star by 
iHeydari-Malaverl ()1988j ). As the optical classification 
is more reliable than our near-IR classification, we will 
adopt a spectral type of 08.5V for #1 as the main ioniz- 
ing star for RCW34. This spectral type is in agreement 
with the Lym an continuum flux derive d from radio con- 
tinuum data (|Caswell fc Havneslll987i h and is also com- 
patible with the He 1 and Br7 emission of the H 11 region 
(Sect. 3.3.1). This corresponds to a mass of 20 M Q for 
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TABLE 3 

Equivalent widths and /C-band spectral types of the OB stars in RCW 34. 



Star 


Brll (1.68 fim) 


He I (1.70 /an) 


He I (2.11 /im) 


Br7 (2.166 /an) 


K-Spectral type 


Optical Spectral type 


1 


2.05 ± 0.3 A 


0.91 ± 0.2 A 


0.4 ± 0.4 A 


4.43 ± 0.4 A 


k09-Bl 


08V - B1V 


2 


2.90 ± 0.4 A 


0.75 ± 0.2 A 


0.6 ± 0.2 A 


6.07 ± 0.3 A 


kB2-B3 


B1V - B2.5V 


3 


5.79 ± 0.7 A 


0.76 ± 0.2 A 


< 0.4 


6.06 ± 0.4 A 


kB4-B7 


B3V - B7V 




Fig. 7. — Observed K vs. (J-K) color-magnitude di agram of th e 
region in RCW 34 covered by SINFONI taken from|firE| 12001 1. 
Stars brigher than K=1A mag (horizontal dotted line) are marked 
with symbols representing their spectral types derived in Sect. 4.1 
and 4.2 ( star symbols: OB stars, diamonds: G stars, triangles: K 
stars). The featureless stars are displayed as crosses. The dashed 
line re presents the 2 Myr isochrone taken from Lcicunc & S chaererl 
l200iT) . 



TABLE 4 

The properties of the detected OB stars in RCW 34 



Star 


Spectral Type 


K (mag) 


J — K (mag) 


A v (mag) 


1 


08.5 V 


8.96 ± 0.01 


0.73 ± 0.01 


4.2 ± 0.1 


2 


B0.5 V - Bl V 


10.35 ± 0.02 


0.68 ± 0.03 


3.9 ± 0.2 


3 


B2 V - B3 V 


11.61 ± 0.03 


0.75 ± 0.05 


4.1 ± 0.4 



star #1 (jMartins et alJl2005f ). 

The measured color of early type stars is a good mea- 
sure of the extin ction. For the intr insic colors we used the 
values give n by iKoornneej (119831) . We used the extinc- 
tion law of lCardelli et al.l (|1989f ); the derived extinction 
towards the three stars is listed in Table [4] and ranges 
from A v =3.9 to 4.2 mag (A K = 0.42 to 0.45 mag). The 
error in the Ay values reflects both the photometric er- 
rors and the uncertainties in ( J-K)o due to uncertainties 
in their spectral classification. These values are in good 
agreement with those derived from the Brl0/Br7 ratio 
of the Hn spectra (Sect. 3.3.1). 

The spectral types and extinctions of the OB stars, 
coupled with their absolute magnitudes, allow a dis- 
tance determination. For th e Q8.5V star #1 , we ob tain 
the absolute mag nitude from iMartins fc Plezl (|2006l ) and 
lAller et al.l (|1982l ). This results in a distance of 2500±200 
pc for the 08.5V star. The uncertainty r eflects the differ- 
ence in the absolute magnitudes between [Martins fc Pled 



(2006) andJAller et al.l (|1982l ). Hanson et all (|1997t ) and 
IBlum et "aTT i 20001 ) provide absolute magnitudes for Zero 
Age Main Sequence (ZAMS) stars. Their results are 
about one magnitude fainter than corresponding values 
for the main sequence of the Hertzsprung-Russell dia- 
gram (HRD). If we adopt the ZAMS magnitudes, we 
obtain a distance of 1800 ± 200 pc. The cluster age is 
probably a few Myrs, hence we prefer the main sequence 
values. Thus, we adopt 2500 ± 200 pc as the distance to 
RCW 34. This distance is slightly closer than the 2.9 kpc 
derived by Heydari-Malayeri (1988) using the same spec- 
tral type and optical magnitudes. The difference results 
from the different absolute magnitudes and extinction we 
apply. 

Using the 2.5 kpc distance, we can refine the spec- 
tral type classification of the two B stars in the cluster. 
Using the distance modulus, the Ay and the observed 
brightness we derive an absolute if-band magnitude of 
M(K) = -2.1 mag for #2 and M(K)= - 0.8 mag for #3. 
With the calibration of lAller et al.l (|1982[ ) this results in a 
B0.5 - Bl V spectral type for #2 and a B2-B3V spectral 
type for #3 (Table 0J, still compatible with the original 
near-infrared spectroscopic classification. These spectral 
types correspond to a mass of 11 and 7 M respectively 
(|Aller et alJfl98l . 

4.2. Late-type stars 

Besides the 3 OB stars we found 21 stars showing ab- 
sorption lines typical of later spectral type (Fig HJ . The 
most prominent lines are the CO first overtone absorp- 
tion bands between 2.29 and 2.45 /im, and absorption 
lines of Mg I and other atomic species in the ii-band as 
well as Cai and Nai in the iv"-band. To determine if 
these stars are cluster members, we need to obtain their 
spectral type and luminosity class. 

For the classification of the late- type st ars w e use 
the reference spec tra of iCushing et al.l ([2005) and 
iRavner et al.l (|2009l ). These spectra are taken with the 
IRTF telescope 16 and cover the spectral types between 
F and M in all luminosity classes with a resolution be- 
tween R = 2000 and 5000. These reference spectra were 
convolved and rebinned to the resolution of SINFONI. 
We corrected for the differences in spectral slope due to 
the extinct ion of our SINFONI t argets using the extinc- 
tion law of lCardelli et al.l (|1989f ). The atomic lines, like 
Mgi and Nai, are used for the temperature determina- 
tion, while the CO and OH lines serve as a luminosity 
indicator. 

Comparison of the SINFONI with reference spectra 
indicates that the CO (2.29 /im) absorption is usually 
deeper than in dwarf reference spectra, but not as deep 
as in the giant spectra. In the few cases where a lu- 
minosity class IV reference spectrum was available, the 

16 http://irtfweb.ifa.hawaii.edu/" spex/WebLibrary/index.html 
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Fig. 8. — Normalized SINFONI H-band (left panel) and _fC-band spectra of the three OB stars detected in RCW 34. Their spectra are 
characterized by absorption lines of the hydrogen Brackett series and He I (dotted vertical lines) . 



spectrum provided a better match to the observed spec- 
trum. This suggests that our late type sta r s are low- 
and intermediate PM S stars; indeed. ILuhmanl ([1999D and 
IWinston et al.l (|2009D find that PMS spectra have a sur- 
face gravity intermediate between giant and dwarf spec- 
tra. If the stars were giants, dust veiling could make the 
if-band CO lines weaker. However, the ii-band CO and 
OH lines are also much weaker, while the atomic lines 
have the expected EWs. Therefore, this seems to be a 
surface gravity effect, suggesting a PMS nature of our 
late type stars. 

We used the relation from lKenvon fc Hartmannl £l995) 
to convert the spectral type into effective temperature. 
This relation applies only for main sequen ce stars; for 
PMS stars a different rel a tion may hold ([Hillenbrand! 
[T997t IWinston et al.ll2009h . ICohen fc Kuhl (119791) show 
that the temperature of PMS stars might be overesti- 
mated by values between 500 K (G stars) and 200 K 
(mid-K). We took this source of error into account when 
calculating the errors in the effective temperature. 

Table [S] shows the result of the classification for all the 
late-type stars brighter than K=14 mag. The spectral 
types range from early G to mid K; about half of the 
stars have a surface gravity indicative of a PMS nature, 
their luminosity class is marked as "PMS". The stars 
classified as dwarfs are candidate foreground stars, but 
this does not exclude them as PMS cluster members, as 
IWinston et al.l (|2009D also find PMS stars with surface 
gravity very close to the dwarf values. 

With the knowledge of the spectral type, the observed 
color can be converted into an extinction estimate. The 
intrinsic J — K color for dwarfs and giants can differ 
as m uch as 0.4 magnitudes for late K stars (jKoornneeJ 
1983). As the late-type stars are PMS stars, we used as 
the intrinsic J — K color the mean of the color for dwarfs 
and giants. The difference between the mean value and 
the dwarf and giant values is used as the error in the 
intrinsic color. The derived extinctions are listed in Table 
[H Source #17 does not have a J-band magnitude as 
the photometry is too contaminated by source #1; it is 
excluded from further analysis. 

All the stars with A v > 12 mag (#6, #8,#16 and 
#24) are located North of the Br7 emission in the re- 
gion where the H2 and IRAC dust emission is present. 



TABLE 5 

T EFF , SPECTRAL TYPE, Av AND IR-EXCESS OF THE 
DETECTED LATE- TYPE STARS BRIGHTER THAN K=li MAG. 



Star 


T e// (K) 


Sp. Type 


Lum. class 


IR excess 


4 


5520 


± 365 


G8 ± 1 


PMS 




5 


5410 


± 365 


G9 ± 1 


PMS 


n 


6 


5080 


± 380 


Kl ± 1 


PMS 


+ 


7 


5860 


± 415 


G2 ± 2 


V 


+ 


8 


5680 


± 360 


G6.5 ± 1 


V 




9 


4900 


± 350 


K2 ± 1 


PMS 


+ 


10 


5520 


± 355 


G8 ± 1 


PMS 




11 


5800 


± 405 


G4 ± 1 


PMS 


+ 


12 


4900 


± 445 


K2 ± 2 


PMS 


n 


13 


5680 


± 360 


G6.5± 1 


PMS 




14 


5680 


± 360 


G6.5 ± 1 


V 


+ 


15 


5520 


± 355 


G8± 1 


V 




16 


5080 


± 380 


Kl ± 1 


PMS 


n 


17 


5080 


± 380 


Kl ±1 


PMS 


n 


IS 


5080 


± 380 


Kl ± 1 


V 


+ 


19 


4900 


± 350 


K2 ± 1 


V 


+ 


20 


4900 


± 350 


K2 ± 1 


PMS 


n 


23 


5520 


± 355 


G8 ± 1 


PMS 


+ 


21 


4800 


± 515 


early K 


? 


n 


25 


4800 


± 515 


early K 


? 


+ 


26 


5410 


± 365 


G9 ± 1 


V 


n 



Note. — '+': infrared excess is present, '-' no infrared 
excess detected, 'n': infrared excess could not be deter- 
mined due to bright background emission. 



Their extinction is consistent with being embedded in the 
molecular cloud North of the H 11 region. The objects sit- 
uated amidst the Br7 emission show a visual extinction 
between 4 and 10 mag, consistent with that measured 
from the nebular line spectrum (Sect. 3.2.1). 

To confirm the PMS nature of the objects we looked 
for Br7 emission from the source. We found two point 
sources in the continuum subtracted Bi"7 line map. One 
of the sources corresponds to #23, a G8 PMS star. Fol- 
lowing [Muierone^OL| (|1998|) . this Br7 emission could 
arise from accreting matter. The other Br7 emitting 
point source is #27 (Table [6]). This star is rather faint 
(K= 15.5 ± 0.1 mag) and the continuum has too low 
S/N to detect any absorption features. The emission 
line spectrum also shows Hei (2.058 ^m) and the Brack- 
ett lines in the ii-band, identical to the surrounding 
H11 region (Fig [5]). The Br7 emission from this ob- 
ject could arise from a disk being photo-evaporated by 
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Fig. 9. — Normalized SINFONI H- band (left panel) and _R"-band (right panel) spectra of the 21 objects showing late-type stellar spectra. 
Most of these objects are PMS stars. The vertical dashed lines show the location of some of the absorption lines used for the classification 
of the stars. The stars are numbered according to Table [6] 
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the main ionizi ng star. Th ese "proplyds" have been 
found in Orion (|O'Dellll2001[) and other high -mass star- 
for ming regions like M8 (|Stecklum et al.lll998ft and RCW 
38 (jDeRose et alJ l2009h . The spatial resolution of our 
data (~0.6", 1500 AU at 2.5 kpc) did not allow us to re- 
solve the proplyd, which is consistent with the si zes (up 
to 50 AU) measured for the Orion Proplyds (jO'Delll 
1998). The observed K-band magnitude and lower-limit 
of source #27 does seem consisten t with near-infra red 
observations of the Orion Proplyds (jRost et al.l 120081 ). 

In the spectra of e.g. source #9 and #17, emission 
and absorption of Br7 is detected as well. These lines are 
artifacts from the Br7 removal using the median filtering. 
In the Br7 line image, no extra emission or absorption is 
present. 

Two objects brighter than if =14 mag are unclassified. 
These objects (#21 and #22) have very red spectra and 
a rather low S/N. We could not identify any features in 
their spectra. As the spectra are extremely red, they are 
likely dominated by dust emission and therefore the spec- 
trum of the underlying star might not be visible. Like 
the late-type stars with high extinction, these feature- 
less spectra are located in the high-extinction area, in 
the North of the SINFONI field and therefore are most 
likely embedded stars. 

The remaining sources are too faint to obtain a proper 
classification. In the spectra of some of these objects, CO 
absorption at 2.29 ^m can still be detected, but the spec- 
tra are of too low S/N to detect the atomic lines needed 
for an effective temperature classification. A handful of 
the objects have a J — K color around 1.0 mag (see Fig. 
[7|), indicating that these might be late- type foreground 
stars with very little extinction. Most of the faint objects 
are concentrated around J — K = 2.0 mag. This color 
is similar to that of the brighter PMS stars around Ay 
= 5 — 10 mag. The redder and more deeply embedded 
objects are not observable as they are too faint to detect 
in the J-band. The objects with K > 14 mag, clustered 
around J — K = 2 mag, are most likely PMS stars with 
masses below 1.5-2 M . 

4.3. IR excess 

In total, 17 sources (out of 26) have been detected in 
one or more Spitzer /IRAC bands with only three sources 
detected in all four bands. Most of the sources are only 
detected at 3.6 and 4.5 ^m. At longer wavelengths, the 
high and variable background emission makes the detec- 
tion of point-sources very difficult. The objects detected 
in all 4 bands are labeled in Fig. [2] One of the PMS 
stars (# 11) is identified as a class II object. The other 
2 objects (# 14 and #21) are classified as class 0/1 and 
will be discussed in Sect. 5.3. 

To determine if the PMS stars have an infrared excess, 
we construct the SEP and com pare it with low-resolution 
Kurucz spectra (jKuruc 3 fl993h for the 1 to 8 y«m range. 
This allows us to determine the IR excess even for ob- 
jects detected in only 1 or 2 IRAC bands. Column 5 of 
Table [5] shows that 9 PMS stars show a detectable IR 
excess (marked as a '+'), while for 5 sources no IR ex- 
cess is detected (marked as a '-'). The remaining sources 
are marked with an " n" : no IRAC counterpart could be 
detected due to blending with other sources or spatial 
coincidence with the bright, extended arc emission. The 
infrared excess is already present at 3.6 /im with an ex- 



cess of more than a magnitude in some cases. The de- 
tected IR excess of at least 9 of the PMS star candidates 
reinforces the hypothesis that these stars are young and 
still surrounded by a circumstellar disk. 

The SINFONI spectra, however, suggest that the near- 
infrared excess of most PMS stars is low, as all these stars 
have absorption lines in their if -band spectra. Should a 
strong IR excess be present these lines would not have 
been detected. The if -band spectrum of #11 suggests 
that some veiling might be present; the CO was still 
stronger than that of the best fitting dwarf reference 
spectrum, but the Nai and Cai lines in the if -band were 
very faint. In the spectra of the two sources where no 
absorption features have been detected (#21 and #22) 
a NIR excess is most likely present. These sources also 
have very red colors and a strong IRAC excess. 

5. STELLAR POPULATION 
5.1. Cluster membership 

To discuss the stellar content in detail, we first need to 
separate the cluster members from the fore- and back- 
ground stars. The 08.5V star (source #1) ionizes the 
H II region and is also responsible for the zone of interac- 
tion with the molecular cloud to the north. Analysis of 
the brightness and spectral type of the B stars (sources 
#2 and #3) shows that they are located at the same 
distance (2.5 kpc) and have a similar extinction as the 
main ionizing source and therefore are considered to be 
cluster members. 

It is less obvious that all the late-type stars are cluster 
members as well. As discussed in Sect. 4.2, most of the 
late-type stars have a surface gravity typical for PMS 
stars. Additionally, about half of the stars have an IR 
excess, supporting the result from the spectroscopy that 
those stars are PMS stars surrounded by an accretion 
disk. The extinction derived for these stars is consistent 
with the extinction derived using the HII region and the 
H2 emission. 

Sources #13 and #15 show no evidence for a PMS na- 
ture nor for an IR excess, but have higher visual extinc- 
tions (5.6 and 7 mag, respectively) than observed from 
the H 11 emission in that area. Their early spectral types 
make it difficult to derive a luminosity class because of 
the CO bands are not present. The high extinction sug- 
gests that they are background stars, but then the stars 
should be (super) giants as dwarfs are too faint to be 
seen as background stars. 

Two more stars have Ay values inconsistent with their 
location: sources #14 and #26 show very low extinction 
(0.25 and 1.8 mag, respectively) while they are located 
towards the high extinction area North of the H 11 region. 
Star #14 has virtually no extinction and a surface grav- 
ity typical for main sequence stars; however, it coincides 
with an IRAC source having colors of a class 0/1 object. 
Also, on the Br7 line map some faint Br7 emission is 
seen surrounding the star. The fact that this star has 
no extinction is inconsistent with its IRAC colors and a 
protostar nature. If the star were a protostar, it would 
be more deeply embedded and have a much larger Ay . It 
is most likely a chance alignment of a foreground late G 
star (as seen in the NIR) with a deeply embedded proto- 
star only seen at mid-IR wavelengths. Additionally, star 
# 26 (Ay = 1.8 mag), also located towards the high- 



14 



Bik ct al. 




08 h 56™ 32 s 30 s 28 s 26 s 24 s 

Right Ascension 



Fig. 10. — Stars brighter than K=14 mag are indicated on the SINFONI 2.15 micron continuum map. Overplotted in gray contours is 
the Br7 line map. The sources are marked with different symbols representing their spectral types. The star symbols represent the 3 OB 
stars, the diamonds the G stars and the triangles the K stars. The featureless stars (#21 and #22) correspond to the cross signs. The 
location of the proplyd (#27) is marked with an asterisk. 




Fig. 11. — Left: The extinction corrected K vs log(T e fr) HRD. The A"-band magnitude has been corrected for the distance modulus. 
The plotting symbols for the stars are the same as in Fig. ]7] The filled symbols represent objects with infrared excess detected with IRAC 
(Table [5}, and the open symbols those without IR excess. Note that in a large fracti on of the objects without I R excess, the background 
emission is very bright. Over-plotted with a do tted line is the 2 My r isochrone from Lcicunc & Schacrer (2001]), and with the solid lines 
the pre-main-se quence evolutionary tracks from Da Rio ct al. (2009). Right: The same data but overplotted with solid lines showing the 
isochrones from Da Rio ct al. (2009). 



extinction area, does not show any PMS characteristic 
nor an IR excess and also is most likely a foreground 
star. 

The locations of the stars in the SINFONI field are 
shown in Fig. 10. The different classes of objects do not 
cover a preferred area of the SINFONI field, but rather 
are spread across the entire field. 

Summarizing, after analysis of their location, redden- 
ing and spectral features we conclude that most of the 



stars are probably cluster members, with two exceptions 
where the stars are probably foreground stars (#14 and 
#26). The nature of sources #13 and #15 is unclear, 
but their high extinction suggests that they are related 
to the cluster. 

5.2. Stellar content and age 

To constrain the mass of the PMS stars and the age of 
RCW 34, we construct a HRD to compare the observed 
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parameters with PMS evolutionary tracks. The derived 
extinctions allow us to de-redden the K-band magnitude, 
convert to absolute magnitude with the derived temper- 
ature, and plot the points in a K vs. log(T G ff) diagram 
(Fig. [ITJ left). We exclude the foreground stars #14 and 
#26 to enable a comparison with the isochrones. The 
symbols reflect the different spectral types. Filled sym- 
bols indicate that a Spitzer/IRAC excess emission was 
detected. 

Over-plotted as a dotted line is the 2 Myr main se - 
quence isochrone taken from lLeieune fc Schaererl l pnOlh . 
The 2 Myr isochrone is chosen to illustrate the main se- 
quence; 1 Myr or 5 Myr isochrones would have served 
equally well. A 20 M star does not evolve signifi- 
cantly off the main sequen ce for the first 8 Myr of its life 
(jLeieune fe Schaererll200lT ). therefore we adopt 8 Myr as 
an upper limit for the age of RCW 34. 

The solid lines in the l eft panel represent the evolu- 
tionary tr acks taken from | Da Rio et al.l (|2009() using the 
models of iSiess et al.l (|2000f) for 4 different masses; 4, 3, 
2, and 1.5 M . Comparison of the location of the PMS 
stars with these evolutionary tracks yields an approxi- 
mate mass varying from ~ 3 M for the brightest stars 
to 1.5 M . 

Over-plotted in the right panel of Fig. [11] are the PMS 
isochrones (solid lines) between 0.1 and 10 Myr. The lo- 
cation of the G- and K stars is consistent with a range of 
ages. The brightest object (#6) lies between the 0.3 and 
0.5 Myr isochrones; however, the majority of the objects 
span the 1-3 Myr isochrones, with the more massive ob- 
jects closer to the 1 Myr isochrone. For the less massive 
objects, the spread in age is larger, most likely due to 
higher uncertainties in the spectral type determination. 
The location of the PMS stars suggests an age of 2 ± 1 
Myr. 

Comparison of the location of the PMS stars in 
the HRD with those of the Herbig AeBe stars 
(|van den Ancker et al.l li~998l ) shows that the PMS stars 
are younger than the Herbig AeBe stars and will evolve 
from their present G- and K spectral type to late B, A or 
early F spectral type when they become main sequence 
stars. This also explains the fact that no A or F stars are 
spectroscopically detected, as they are still in the PMS 
phase at cooler temperatures. 

Even though the PMS stars are younger than the Her- 
big AeBe stars, they h ave, in contrast to the Herbig stars 
(e.g. lEiroa et al.ll2001l ). no K-band excess. This could be 
the result of a different disk structure, but also caused 
by the fact that the PMS stars are emitting most of their 
energy in the near-infrared wavelength regime, where the 
SEDs of Herbig AeBe stars peak in the optical, thereby 
outshining the contribution of the circumstellar disk. 

The extinction measurements show that the sources 
North of the H n region are more deeply embedded in 
the molecular cloud. A younger age for this region is 
suggested by the presence of several protostars as well as 
outflows and a dense molecular cloud, which would sug- 
gest a younger age for the PMS object in this area as well. 
The PMS stars North of the Hll region (# 6, #8, #16, 
#24) all have the brightest absolute magnitude in their 
temperature bin with source # 6 as an extreme example. 
Their position in the HRD could suggest that they are 
the youngest of the stars of given spectral type. However, 
their high-extinction leads to a selection effect (only the 



brightest sources are detected). Fainter sources would 
fall below the limit where SINFONI spectra could be 
classified. So, if there were an age difference, the current 
data would not be able to quantify this. Deeper spec- 
troscopic observations to find more of the high-extintion 
PMS stars North of the H II region would be required to 
make any definitive statements. 

5.3. Deeply embedded population 

As the SINFONI observations show, large extinction 
variations are present in RCW 34. The high extinction 
towards the North means that we miss sources that are 
too deeply embedded to be visible in the K-band. In 
our IRAC catalogue, we found 2 sources with no near- 
infrared counterpart (#28 and #29 in Table [6]). Both 
sources have IRAC colors of deeply embedded protostars 
(class 0/1). 

Source #28 is located in the North-East corner of the 
SINFONI field, only 3" away from one of the outflows 
detected in the H2 line emission (outflow c, Sect. 3.2.2). 
The other embedded source (# 29) is located very clo se to 
source #1, but is not related. I Thomas et al.l ()1973l ) and 
iWilliams et al.l (|1977l ) detected #29 as a mid-infrared ex- 
cess towards vdBH 25a. The IRAC images show, how- 
ever, that it is a separate source. The SED of # 29 is 
difficult to determine; at 3.6 /im the contribution of the 
08.5V star is still present, and at 5.8 and 8.0 /im this 
object is the brightest source at MIR wavelengths. To 
determine a better SED, high- resolution and deeper mid- 
infrared observations are needed to discriminate between 
the emission from the protostar and the arc in the north. 
Additionally, as discussed above, the IRAC counterpart 
of source #14 could be another embedded protostar not 
seen at NIR wavelengths. 

Outside the SINFONI field, several other protostars 
are located. Some of these objects seem to ha ve a K-b&nd 
counterpart in the SOFI image oflBiia ([200l . suggesting 
that the objects in our SINFONI field are even more 
embedded than those located a bit further away from 
the arc. 

5.4. Stellar population of the surroundings 

IRAC imaging reveals a large bubble surrounded by 
extended 8 /im emission. Inside this bubble is a cluster 
of intermediate and low-mass class II objects. Typical 
ages for class II objects are of the order of a few Myrs 
(|Kenvon fc Hartmann|[T995H White et al.ll2007ft . not dis- 
tinguishable from the age of the stars in the H 11 region. 

The class 0/1 stars, identified mainly close to the bright 
arc at the northern edge of the Hl l regio n, might be sig- 
nificantly younger. iWhitnev et al.l ((2003) sugg est an age 
of ~ 1 5 years for class I sources; however, IWhite et al.l 
(|2007t) present evidence that the low-mass class I objects 
in Taurus have a similar age as the class II sources. The 
location of the class 0/1 sources coincides with methanol 
masers, CS emission, high extinction and outflows, and 
therefore are most likely younger than the objects inside 
the bubble and the H 11 region. 

Comparing the 2MASS if-band magnitudes of the 
class II stars in the bubble with the if-band magnitudes 
of the PMS stars in the Hn region shows that 13 out of 
17 class II objects are about a magnitude brighter. The 
other four are of similar magnitude. Additionally, in a 
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(J — H) vs. (H — K) color-color plot, most of the class II 
sources have a near-infrared excess and cover the same 
region in the color -color diagram as the Herbig Ae stars 
(jEiroa et al.ll200l . 

The difference in if-band magnitude could be ex- 
plained with the presence of a -ftT-band excess, suggesting 
that the class II sources are more similar to the Herbig 
Ae stars and are therefore more evolved than the PMS 
stars. An optical spectroscopic classification is needed to 
confirm this. 

The location of the central star of the HII region, 
near the northern edge of the bubble suggests that this 
source might not have cre ated the bubble. Also the Ha 
(|Parker fc Phillippsl I1998D as well as the Br 7 emission 
peaks at the location of the O star. In the rest of the 
bubble only faint Ha emission is present, indicating the 
absence of massive stars inside the bubble. Apart from 
the class II sources detected in the bubble, one source in 
the center (Fig. [3]) is detected without infrared excess. 
The magnitude at 3.6 fim (10.0 ± 0.02 mag) is compat- 
ible with a BOV star at the distance of RCW 34. How- 
ever, the optical p hotometry of the source taken from the 
USNO catalogue (|Monetlll998f ) , is not compatible with a 
(reddened) early type star. Therefore we conclude that 
this star an unrelated, foreground object. 

This raises the question whether other mechanisms 
could be responsible for the bubble. The scenario that 
the bubble would be created by a supernova seems un- 
likely. As the stars in the bubble are only a few Myrs old, 
the supernova needs to come from an early type O star. 
Assuming a normal IMF, the other O stars still would 
have to be visible. 

Massive stars affect their environment by ionizing pho- 
tons and stellar winds. Estimating the mechanical lumi- 
nosity (E mec /j = 1/2 Mv 2 ) needed to blow a bubble with 
a radius of 2.5 pc inside a medium with a typical molec- 
ular cloud density (10 3 cm -3 ), expanding at 10 km/s 
would be on the order of 10 42 erg/s. As a comparison, 
the wind lu minosity of a mid-B star is around 10 30 erg/s 
(|Tanll2004n . 

Assuming ^10 44 ionizing photons s , typical for a 
mid-B star and a density of 10 3 cm -3 , the Stromgren 
radius is 0.014 pc. This is a factor 100 lower than the 
observed radius of the bubble (2.5 pc). Assuming that 
the 20 observed class II stars all emit the same number 
of ionizing photons and assuming a density of 100 cm" 3 , 
the Stromgren radius would be comparable to that of the 
measured bubble size. Additionally, the Hn will expand 
into the molecular material due to the higher pressure, 
making this the most plausible scenario. 

Another possibility is that the OB star which created 
the bubble has been expelled from the center of the clus- 
ter due to a binary interaction, similar to what is sug- 
geste d for the origin of BN in Orion ((Rodriguez et al.l 
2005) . The central star of the H n region could be one of 
the stars from the disrupted binary, with the other star 
moving in the opposite direction. With the current data 
we cannot exclude this scenario, detailed velocity mea- 
surements of the star in the HII region could give more 
insights into this scenario. 

6. DISCUSSION 

The interaction region between the H n region and the 
molecular cloud North of it raise the suggestion that 



some form of triggered star formation plays a role here. 
After having derived several properties of the H II region, 
the molecular cloud and the adjacent bubble we can now 
attempt to reconstruct the way this region has formed. 

RCW 34 is composed of a large bubble in which 
intermediate- and low-mass PMS stars have formed with 
no evidence for a central ionizing source. The ionizing 
source and its surrounding Hn region, however, are lo- 
cated at the edge of the bubble, very close to the molec- 
ular cloud. The H n region and the molecular cloud are 
separated by an ionization front and a PDR traced by 
H2 emission and 5.8 and 8.0 /xm emission. The age of 
the stars inside the bubble (several Myr) is not known 
with sufficient accuracy to detect an age difference with 
the stars inside the H 11 region (2 ± 1 Myrs) . The age of 
the objects forming inside the molecular cloud, however, 
is most likely younger, as several tracers of active star 
formation have been detected. 

Inside the H 11 region, a North-South velocity gradient 
has been detected from Ha spectroscopy, res embling that 
of S88B, a proto- typical c hampagne flow (|Garav et al.l 
[19981: IStahler fc Pallal[2005h . The detection of the cham- 
pagne flow, as well as the presence of the ionization front 
and PDR in the North, suggest that the Hn region is 
formed at the edge of the molecular cloud. Due to the 
lower-density in the bubble, the HII region can expand 
into the bubble creating the observed velocity gradient. 

We already can exclude one of the triggered star- 
formation scenarios. The collect and collapse scenario, 
where the star-forming clumps are formed by swept-up 
material is very likely not the scenario that plays a role 
here. Firstly, only at the northern edge of the bubble is a 
molecular clump pres ent, while more clumps around the 
bubble are expected (jDeharveng et al.l 120051 ). Secondly, 
it is very unlikely that the energy output of a cluster of 
only intermediate mass stars would be enough to sweep 
up a dense molecular clump of more than 1000 M Q as 
detected in CS emission. 

Also unlikely is that the formation of the Hn region 
and the O star caused both the bubble and the star for- 
mation in the Molecular Cloud. The difference in evolu- 
tionary stage (class II stars in the bubble and class 0/1 
stars in the molecular cloud) is difficult to explain by the 
triggering due to the H 11 region alone. 

The location of the bubble, Hn region and molecular 
cloud suggest a sequence of events. First, the cluster 
of stars formed which created the bubble. After that, 
the HII region formed at the border of the bubble and 
the larger molecular cloud; the latter currently forming 
stars. The location of the Hn region with respect to 
the molecular cloud as well as the presence of the PDR 
shows that the O star is influencing the molecular cloud 
and probably also inducing the current star formation. 

The question remains whether the formation of the 
bubble triggered the other star-formation events. A 
mechanism whereby the energy output of the stars inside 
the bubble would be pushing a pre-existing core to col- 
lapse could be at work here. If a molecular core present 
at the edge of the molecular cloud were in equilibrium, 
a small push from the bubble could have made the core 
collapse and enable the formation of the OB stars. 

Another alternative would be that the triggering mech- 
anism is external and that a wave of star formation is 
propagating through the molecular cloud resulting in the 
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formation of the different sites of star formation. 

To determine which scenario is the correct one, a better 
age determination of the class II sources by means of 
optical or near-infrared spectroscopy would be needed. 
A triggering scenario would result in an age difference 
between the stars in the bubble and the H II region. 

7. CONCLUSIONS 

In this paper we have presented SINFONI integral 
field spectroscopy of the H n region RCW 34 as well as 
Spitzer /IRAC imaging of the surrounding area. The re- 
sults can be summarized as follows. 

1. RCW 34 consists of 3 separate regions, a large bub- 
ble outlined by Spitzer /IRAC emission with an Hii re- 
gion at the northern edge of the bubble and a dense 
molecular cloud just North of the H n region. IRAC pho- 
tometry of the sources inside the bubble reveals a cluster 
of intermediate class II sources and shows that no OB 
star is located near the center of the bubble. The class 
0/1 objects are mostly detected towards the dense molec- 
ular cloud in the north, suggesting active star formation 
present in the molecular cloud. 

2. SINFONI spectroscopy of the Hii region reveals 
an extinction gradient from North (Ay =15.9 mag) to 
South (Ay = 1.5 mag), supportive of the champagne 
flow hypothesis measured in the Ha velocity profile. The 
region between the Hii region and the molecular cloud 
is identified as a PDR, where the H 2 emission originates 
from UV fluorescence. Additionally, two outflows are 
identified based on their H2 emission. 

3. SINFONI spectroscopy identifies the 3 brightest 
stars as ionizing sources of the Hn region. The main 
ionizing source, #1, is of spectral type 08.5V while the 
two other stars are B0.5V and B3V stars. This results in 
a revised distance estimate of 2.5 ± 0.2 kpc for RCW34. 
Additionally, we detected 21 G and K stars of which two 
stars are foreground objects and 19 objects are PMS clus- 
ter members. With a mass between 2 and 3 M , these 
objects are still contracting to the main sequence and are 
the precursors of the A and F stars. Combining the lo- 
cation of those objects in a HRD with PMS evolutionary 
models results in an age estimate of 2 ± 1 Myr for RCW 
34. 

4. In this scenario, the cluster stars created the bub- 



ble, which then collapsed a pre-existing dense core at the 
edge of the molecular cloud. The resulting H 11 region 
would encounter the density gradient at the edge of the 
cloud, and thus develop as a champagne flow. Whether 
the bubble was the original trigger, or if there was some 
other, external trigger, cannot be answered with the cur- 
rent data; a detailed spectroscopic study of the bubble 
stars will be required. 

This study shows that the interaction between the 
young massive stars and the surrounding molecular cloud 
could not only lead to the destruction of the cloud, but 
also to the formation of a younger generation of stars. 
In the majority of our regions we have identified mul- 
tiple generations of star formation. Several papers are 
in preparation to discuss the remainder of the sample of 
10 clusters observed with SINFONI (e.g. Wang et al, 
in prep, Vasyunina et al. in prep.) and more final con- 
clusions will be presented as a final paper discussing the 
survey. 
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TABLE 6 

Photomeric and spectroscopic properties of the sources located in the SINFONI field of view. 



OLcLI 










[3.6J 


|4.t>J 


[b.aj 


[8.0] 




A Tr 


Comments 




(fl ni s) 


jo 1 ll\ 

\ ) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(mag) 




1 


08:56:28.1 


-43:05:55.6 


8.96 ± 0.01 


0.73 ± 0.01 










08.5V 


4.8 ± 0.1 


vdBH 25a 


2 


08:56:28.1 


-43:05:57.1 


10.35 ± 0.02 


0.69 ± 0.03 










B0.5-B1V 


4.1 ± 0.2 




3 


08:56:29.8 


-43:05:59.4 


11.61 ± 0.03 


0.76 ± 0.05 


10.81 ± 0.03 


10.37 ± 0.06 






B2-B3V 


4.3 ± 0.4 




4 


08:56:31.7 


-43:05:56.3 


11.66 ± 0.04 


1.56 ± 0.06 


11.26 ± 0.03 


11.01 ± 0.06 






G8±l PMS 


4.7 ± 0.6 




5 


08:56:28.2 


-43:05:58.3 


11.73 ± 0.04 


1.54 ± 0.06 










G9±l PMS 


4.6 ± 0.6 




6 


08:56:29.1 


-43:05:31.6 


12.01 ± 0.04 


3.53 ± 0.15 


9.67 ± 0.02 




8.51 ± 0.16 




Kl±l PMS 


14.2 ± 0.9 




7 


08:56:29.6 


-43:05:51.9 


12.23 ± 0.05 


1.42 ± 0.07 


10.84 ± 0.17 








G2± 2 V 


4.4 ± 0.6 




8 


08:56:31.0 


-43:05:32.4 


12.61 ± 0.05 


3.6 ± 0.2 


11.61 ± 0.07 








G6.5 ± 1 V 


15.0 ± 1.3 




9 


08:56:29.7 


-43:06:03.2 


12.68 ± 0.06 


1.93 ± 0.11 


10.88 ± 0.03 


10.49 ± 0.04 


9.7 ± 0.10 




K2 ±1 PMS/V 


5.7 ± 0.7 




10 


08:56:29.7 


-43:06:24.3 


12.83 ± 0.06 


1.72 ± 0.11 


12.17 ± 0.03 


11.88 ± 0.05 






G8 ±1 PMS 


5.5 ± 0.8 




11 


08:56:31.9 


-43:06:01.5 


12.87 ± 0.06 


2.71 ± 0.16 


11.00 ± 0.02 


10.36 ± 0.02 


9.76 ± 0.04 


8.92 ± 0.08 


G4 ±1 PMS/V 


10.8 ± 1.1 


class II 


12 


08:56:29.1 


-43:06:19.1 


12.88 ± 0.06 


1.80 ± 0.11 










K2 ±2 PMS/V 


5.0 ± 0.7 




13 


08:56:27.8 


-43:06:21.4 


12.99 ± 0.06 


1.99 ± 0.13 


12.17 ± 0.04 


11.67 ± 0.07 






G6.5 ± 1 V 


7.0 ± 0.9 




14 


08:56:26.8 


-43:05:35.2 


13.01 ± 0.07 


0.71 ± 0.09 


10.51 ± 0.12 


9.89 ± 0.10 


7.10 ± 0.11 


5.21 ± 0.11 


G6.5 ± 1 V 


0.3 ± 0.7 


vdBH 25a C 


15 


08:56:26.7 


-43:06:27.0 


13.09 ± 0.07 


1.73 ± 0.12 


12.38 ± 0.04 


12.00 ± 0.07 






G8 ± 1 V 


5.6 ± 0.9 




16 


08:56:29.2 


-43:05:32.5 


13.22 ± 0.07 


3.32 ± 0.25 










Kl ± 1 PMS 


13.1 ± 1.5 




17 


08:56:28.0 


-43:05:55.2 


13.23 ± 0.07 


• • ■ ± 0.08 










Kl± 1 PMS/V 






18 


08:56:29.4 


-43:06:00.3 


13.29 ± 0.07 


2.17 ± 0.16 


10.90 ± 0.03 


10.21 ± 0.04 


9.42 ± 0.10 




Kl ± 1 V 


7.2 ± 1.0 




19 


08:56:27.8 


-43:06:07.5 


13.47 ± 0.08 


1.60 ± 0.14 


12.01 ± 0.10 


11.06 ± 0.11 






K2 ± 1 V 


4.0 ± 0.9 




20 


08:56:24.6 


-43:06:01.1 


13.54 ± 0.08 


2.66 ± 0.21 










K2 ± 1 PMS 


9.4 ± 1.3 




21 


08:56:29.6 


-43:05:34.1 


13.70 ± 0.09 


> 4.20 


10.01 ± 0.03 


9.03 ± 0.03 


8.23 ± 0.08 


6.96 ± 0.13 


featureless 


class 0/1 




22 


08:56:30.7 


-43:05:35.4 


13.76 ± 0.09 


> 4.14 


11.25 ± 0.07 


10.53 ± 0.06 






featureless 




23 


08:56:28.0 


-43:06:13.2 


13.76 ± 0.09 


2.05 ± 0.19 


12.40 ± 0.11 


11.66 ± 0.15 






G8 ± 1 PMS/V 


7.2 ± 1.2 


Br7 emission 


24 


08:56:28.4 


-43:05:34.9 


13.80 ± 0.09 


4.01 ± 0.46 










early K 


15.9 ± 2.8 




25 


08:56:29.1 


-43:06:18.5 


13.88 ± 0.10 


1.57 ± 0.17 


11.75 ± 0.04 


11.27 ± 0.06 


10.8 ± 0.15 




early K 


3.6 ± 1.4 




26 


08:56:26.8 


-43:05:40.8 


13.93 ± 0.10 


1.06 ± 0.15 










G9 ± 1 V 


1.8 ± 1.0 


vdBH 25q 


27 


08:56:29.0 


-43:05:54.0 


15.5 ± 0.1 


>2.3 














Proplyd 


28 


08:56:32.5 


-43:05:34.69 






11.69 ± 0.03 


9.77 ± 0.02 


8.16 ± 0.02 


6.00 ± 0.03 






class 0/1 


29 


08:56:28.0 


-43:05:53.76 






8.24 ± 0.01 


7.57 ± 0.02 


6.43 ± 0.03 


4.49 ± 0.04 






class 0/1 



C/2 

O 

B 



c 
B 

B' 

4^ 



